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The integrals in Gradshteyn and Ryzhik. 
Part 14: An elementary evaluation of entry 3.411.5 

Tewodros Amdeberhan and Victor H. Moll 



Abstract. An elementary proof of an entry in the table of integrals by Grad- 
shteyn and Rhyzik s presented. 



1. Introduction 

The compilation by I. S. Gradshteyn and I. M. Ryzhik [6) contains about 600 
pages of definite integrals. Some of them are quite elementary; for instance, 4.291.1 

(1.1) / V _ ' dx = — 



x 12 

is obtained by expanding the integrand as a power series and using the value 

' (-I)*" 1 _ TT 2 



(1-2) E 



k 2 12' 

fc=i 

The latter is reminiscent of the series 

fc=l 

The reader will find in [3] many proofs of the classical evaluation (|1.3[) , 

Most entries in [6] appear quite formidable and their evaluation requires a variety 
of methods and ingenuity. Entry 4.229.7 

(1.4) / lnlntanxdx = -ln -^V27r 

Jtt/4 2 \r(i) 

illustrates this point. Vardi |13j describes a good amount of mathematics involved 
in evaluating (jl.4p . The integral is first interpreted in terms of the derivative of the 
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L-function 

MS 

TFT 

(1.6) / In In tan x dx = — j- + L'(l). 

Here 7 is Euler's constant. Then 1/(1) is computed in terms of the gamma function. 
This is an unexpected procedure. 

Any treatise such as [6] , containing large amount of information is bound to have 
some errors. Some of them are easy to fix. For instance, formula 3.511.8 in the sixth 
edition [5] reads 

(17) r dx M)fc 

Jo cosh 2 x ^A/SfcTT' 

The source given for this integral is formula BI(98)(25) from the table by Bierens de 
Haan [2], where it appears as 



The change of variable t — y/x yields a correct version of (11.71) : 

(1.9) / — rr^ = ^Y 



cosh(x 2 ) V2fc + 1 



/ -T2- = L 



It is now clear what happened: a typo produced (|1.7[) . In the latest edition of the 
integral table [6] , the editors have decided to replace this entry with 

dx 
cosh 2 

The right-hand side of (|1.7[) has been corrected. 

The advent of computer algebra packages has not made these tables obsolete. The 
latest version of Mathematica evaluates (jl.lOj) directly, but it is unable to produce 
(Ol) . 

Most of the errors in [6] are of the type: some parameter has been mistyped, 
an exponent has been misplaced, parameters are mistaken to be identical (a common 
mishap is \x and u appearing in the same formula). Despite of this fact, it is a remark- 
able achievement for such an endevour. The accuracy of [6] comes from the effort of 
many generations, beginning with [9] and also including [71 111] . 

A different type of error was found by one of the authors. It turns out that entry 
3.248.5 of 6 is incorrect. To explain the reason for looking at any specific entry 
requires some background. The stated entry 3.248.5 involves the rational function 
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and the result says 
(1.12) 



dx 



(1 + ^)3/2 \ip{x) + yfp{x) 

The encounter begins with the evaluation of 
(1.13) N QA (a;m) = 

in the form 



V 2 2 V6' 



dx 



(1.14) 



N QA (a; m) = - 



(x 4 + 2ax 2 - 
P m (a) 



Um+1 



2 [2(a+l)] m +V2' 

where P m (a) is a polynomial of degree m. The reader will hnd in [H HOj details about 
(|1.13[) and properties of the coefhcients of P m . It is rather interesting that No^a, m) 
appears in the expansion of the double square root function 

- ( _ 1)fc -i 



(1.15) 



, + yTTc= v / oTT + 



l 



fc=i 



-A , 4 (a;fc- l)c fc 



Browsing [5] on a leisure day, and with double square roots in our mind, formula 
3.248.5 caught our attention. After many failed attempts at the proof, a simple nu- 
merical integration showed that (|1.12p is incorrect. In spite of our inability to evaluate 
this integral, we have produced many equivalent versions. The reader is invited to ver- 
ify that, if a(x,p) :— y/ x 4 + 2px 2 + 1 then the integral in (jl.121) is 7(|, 1), where 

dx 



(1.16) 



I(p,q) 



h V a(x, p) y/ a(x, q) y/ a{x, p) + <t(x, q) 
The correct value of \1.\2\ has eluded us. 

The reader is surely aware that often typos or errors could have profound conse- 
quences. In a letter to Larry Glasser, regarding (I1.12[) . we mistyped the function ip(x) 
of ([i~TT|) as 

(1.17) , w = 1 + _i£l_. 

Larry, a consumated integrator, replied with y/3 ^Tanh _1 \/2a; — ^Tanh -1 y/u^j where 

uj = (y/7 - y/3) /2y/7. Beautiful, but it does not help with (fL1"2"]l . The editors of [5] 
have found an alternative to this quandry: the latest edition [6] has no entry 3.248.5. 

Another example of errors in [6] has been discussed in the American Mathematical 
Journal by E. Talvila [12] . Several entries, starting with 3.851.1 [5] 



(1.18) 



f°° b !~tF 

I x sm(ax 2 ) sin(2bx) dx = — \ — 
Jo 2a V 2a 



are shown to be incorrect. This time, the errors are more dramatic: the integrals are 
divergent. These entries do not appear in the latest edition [6]. 
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The website http: //www. math, tulane . edu/^vhm/Table . html has the goal to 
provide proofs and context to the entries in |6|. The example chosen for the present 
article is taken from Section 3.411 consisting of 32 entries. The integrands are combi- 
nations of rational functions of powers and exponentials and the domain of integration 
is the whole real line or the half line (0, oo). There is a single exception: entry 3.411.5 
states that 

f ln2 xdx 7T 2 

(L19) 1 i^ = s 

The next section presents an elementary proof of (|1.19p . 

2. A reduction argument 

The expansion of the integrand in (11.191) as a geometric series yields 

oo 

(2.1) - - =x + Yxe~ kx . 

k=l 

Term-by-term integration produces the following expressions 

/■a j i 00 -ak 00 i 00 p -ak 

(2.2) / -J^-dx = ±a 2 -Y e — +y ^- a y e — . 

y ' In 1 - e~ x 2 ^ k 2 ^ k 2 ^ k 

Jl) fc=l k=l fc=l 

The complexity of these three series decreases as one moves from left to right. We 
now compute each term in (12.2[) . individually. 



The third series. Integrating the geometric series x k = 1/(1 — x) yields ^2 = 



fc=0 n>l 

ln(l — x), which is valid for |x| < 1. Evaluating at x = e~ a gives 

oo p-ak 

(2.3) — = Ml - ^ a ). 

fe=i 

The second series. The Riemann zeta function 
( 2 - 4 ) C(«) = E A 

n=l 

plays a prominent role in the evaluation of the remaining 31 entries in Section 3.411. 
Indeed, the first of these 

is the classical integral representation for (,{v). It is becoming that the special value 

(2-6) C(2) = y 

appears as the second series in 
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The first series. The second series in (|2.2I) is the only remaining part, we are alluding 
to the function x k /k 2 evaluated at x = e~ a . This is the famous polylogarithm 

studied by Euler. See the introduction to [8| for a historical perspective. Aside from 
the series representation 

oo 

(2.7) PolyLog(2,x)=^^, 

k=l 

there is a natural integral expression 

(2.8) PolyLog(2, x) = - f ^ - ^ dt. 



Therefore, the identity (|2.2I) reduces to 

(2.9) £ = \a 2 - PolyLog[2, e~ a ] + ^ - a ln(l - e""), 
and entry 3.411.5 corresponds to the special value 

(2.10) p o iy L0 g[2,I] = ^-iln 2 2. 
Equivalently, 

(2.11) y_J_ = ZL_ m 2 2 . 

v > 2 fc -!fc 2 6 

fc=i 

Euler proved the functional equation 

(2.12) PolyLog[2, x] + PolyLog[2, l-x] = — - lna:ln(l - x) 
for the polylogarithm function. In particular, the case x = h gives (|2.11j) . 

3. An elementary computation of the first series 

A series for In 2 2 can be obtained by squaring In 2 = — so that 



In 2 2 



I n2 n J X I ^ m2 m j ^ nm2"+ r 

\n— 1 / \m— 1 / n.m^l 



v / i - 1 | \ j 



EE 



' (r — m)m / 2 r 

The partial fraction decomposition — = — ( 1 I gives 

(r — m)m r \m r — m ' 

oo 

(3-D ln22 -E^ 

r=l 
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where H r -i = 1 + t; + ■ ■ ■ + — is the harmonic number. Therefore, 

OO ^ OO jj OC 

In 2 + fc 2 — r2 r_1 2 fe_1 fc 2 

fc=l r=l k=l 



r2 r 

r=l 



oo „ 

r2 r - 1 ' 

It remains to verify that this last series is in fact C(2). 

The representation of the harmonic number as 

f 1 1 — x r 

(3.2) H r = / dx 

Jo 1 - x 

gives the desired step. Indeed, if c r is a sequence of real numbers and a is fixed, then 

v /"^ 1 X 

(3.3) ~S~^c r H r a r = I ~ x r )c r a r dx. 

r=l JU r=l 

Thus the function 

OO 

(3.4) f(x)=Y,CrX r 

r=l 

appears in the integral representation 

(3.5) f>iw= r f{a \- f(aX) dx. 



In the present case, c r = l/r2 r 1 and /(a:) = — 21n(l — x/2). Therefore, 

(36) ^__j; ^ +l H r m ^ 2 j^ dv . 

This last integral is computable via (ll.ip and we have come full circle. 

The technique described above in exhibiting an elementary proof of (|2.1ip can be 
applied to 

OO 

(3.7) E^^CO). 

n=l 

J. Borwein and D. Bradley [4] have given 32 proofs of this charming identity. 
Acknowledgements. The second author acknowledges the partial support of nsf-dms 0713836 
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